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Abstract 

We present a semiclassical theory of spin-diffusion in a ferromagnetic metal subject to a temperature gradient. Spin-flip 
scattering can generate pure thermal spin currents by short-circuiting spin channels while suppressing spin accumulations. 
A thermally induced spin density is locally generated when the energy dependence of the density of states is spin polarized. 
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1. Introduction 

Thermoelectric properties attracted new interest in re- 
cent years due to the improved performance of nanometer- 
scale structures [TJ El El 0]. Spin caloritronics addresses 
the interplay of conduction electron charge, spin and en- 
tropy/energy transport in solid state nanostructures in the 
search for improved heat control. The thermodynamics 
of spin flow in metallic structures involving magnetic ele- 
ments is an important research area. The giant magneto- 
thermoelectric power in multilayered nanopillars [5 , ther- 
mally excited spin-currents in metals with embedded fer- 
romagnetic clusters [B], thermal spin-transfer torques [7] 
and magneto-Peltier cooling [5] in spin-valve devices are 
examples of spin-dependent thermoelectric phenomena on 
a nanometer scale. Our understanding of the coupling 
between spin and thermoelectric transport is still very in- 
complete, however. 

Power dissipation is a major problem obstructing further 
miniaturization of electronics, thus deserves to be studied 
in magnetic systems. The Seebeck effect refers to con- 
version of a thermal differential into an electric voltage, 
for instance in a thermocouple. A spin analogue to the 
Seebeck effect has been discovered by Uchida et al. [3] 
in permalloy thin films subject to a temperature gradient. 
The inverse spin Hall voltage picked up by platinum top 
contacts is evidence of a thermally induced spin current 
locally generated by the ferromagnet, which changes sign 
in the middle of the wire and has an approximately linear 
profile between the heat source and drain. 

Here we discuss temperature-gradient-driven thermo- 
electric transport in disordered ferromagnetic metals. 
Starting from the Boltzmann equation we study spin and 
heat diffusion through the ferromagnet in the presence 
of spin-conserving and spin-flip impurity scattering. We 
show that thermal generation of a spin accumulation is 
limited by the spin-flip scattering. At nonzero (finite) 



temperatures the spin-flip scattering equilibrates the spin 
imbalance in the electron distribution function (or spin 
accumulation), but not the temperature-dependent mag- 
netization density at local thermal equilibrium. 

We start in Sec. 2 by formulating a diffusion theory for a 
ferromagnetic metal in the presence of spin-dependent ex- 
ternal forces and derive spin diffusion equations for the dis- 
tribution functions. Spin dependent thermoelectric trans- 
port is studied in Sec. 3 for a ferromagnetic metal subject 
to a temperature gradient. In Sects. 4 and 5 we explain 
the effect of a temperature dependence of the chemical po- 
tential. The spin densities and exchange splitting in the 
presence of temperature-dependent densities of states are 
discussed in the Appendix. 

2. Spin Boltzmann and diffusion equations 

We consider a homogeneous ferromagnetic metal in a 
simple Stoner model with energies ej^ = + aA — fiQ for 
spins a — ±1, where is the kinetic energy for momen- 
tum k — |k|, A (T) the temperature dependent exchange 
splitting, and /i the ground state chemical potential. The 
Boltzmann equation describes evolution of the distribution 
of conduction electrons in the presence of external fields 
and random scattering processes and for a steady state 
reads 



f( q ) 



V k )/ (Q) (k,r) 



(1) 



where = Vk£fc is the electron velocity and /( Q '(k,r) 
its distribution as a function of momentum k, position r 
and spin a. The forces F( Q ) (r) = eV r 0- a (dA/dT) V r T 
originate from the gradients of the electric potential (j) 
and exchange splitting A, where the latter may include 
an external (Zeeman) magnetic field, e — \e\ is the mod- 
ulus of the electron charge and h — h/2n Planck's con- 
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stant. The collision integral on the right-hand side con- 
tains spin-conserving and spin-flip impurity scattering. In- 
elastic scattering is assumed to be strong. The isotropic 
parts of the distribution functions in momentum space 
are then locally thermalized to the Fermi-Dirac form at 
temperature T (r) and spin-dependent chemical potential 
shifts ^ (r) : 



/W( £fe! r) = (47r)- 



dn k /< a >(k,r) 



-M (a) (r))/(fc B T(r)) 



(2) 
(3) 



Inelastic scattering processes then do not appear explicitly 
in the collision term of the right-hand side of Eq. ([!]), 
which in the relaxation time approximation reads (note 
the difference with [TP]): 



dt 



g (a) (k,r) 
fW(e k ,r) -/("«) (e k + 2aA,r) 



.(«) 



(4) 



where ff < a >(k,r) = f<- a Hk,r) - J^\e k ,v) and l/r t ( r a) ( £ ) 



1/tq (e) + 1/t f a> (e) are the relaxation rates that include 
both the spin-conserving and spin-flip scattering times Tq 
and Tj 1 ', respectively. 



(a) 



(«) 



Using fi _1 Vk/^ = v^d^ft") and averaging Eq. (JIJ 
over momentum directions we find for the anisotropic part 
of the electron distribution function 



5 W (k,r) 



(a) (a) 



Substituting Eq. ^ into the Boltzmann Eq. Q, inte- 
grating over k-directions yields the spectral equation 



V 



(a) 



(a) _ 



f (a) _ /(^ Q ) 



.(«) 



where 2?( a ' = v\t^ /Z is the diffusion constant. Intro- 
ducing the spectral current density 



:(«) 



1 



(e,r) = -aM( £ )(v 



F^| £ ) (7) 



the charge current 
J(«)(r)=- e | ^^)/W( k; r) = J<k$ a \e,T) (8) 
and the energy current 



r(") 



(r) 



(9) 



are expressed in terms of the conductivity cr*-"* 1 = 
e 2 7V (Q) X» (a) , in which W (q) is the density of states. 



The charge and energy currents must be conserved, i.e. 
£ (a) V r • jW(r) = and £ (q) V r • j£>(r) = 0. 

The problem is simplified by linearizing Eq. (j6j. We 
disregard the higher order term J?( a )g( a ) and the temper- 
ature dependence of the material parameters. Expanding 
/( a ) around the equilibrium distribution function /g Q \ we 
obtain for the local distribution function in the spin-flip 
terms 

/ (Q) ( £fe ,r)| T =co ns t. « /^(e fe ) + (r) , 



3£i 



and for the divergence 



V r / (Q) (£ fc ,r) 



(10) 




<9A\ 



(a) 



(11) 



After integration over momenta and invoking the Sommer- 
feld expansion, we obtain the particle diffusion equation 



s5 (a) T 



(a) 



(-«) 



(12) 



where l { s f = \Jv^)t^\ ^ = ^ ~ e<j>, and the spin- 
dependent Seebeck coefficients = — eCoTd e In cr^l^ 
(Mott's formula) in terms of the Lorenz constant £ = 
(7r 2 /3) (fcs/e) 2 have been introduced. By multiplying Eq. 

Wat by e k and integrating, we obtain the heat diffusion 
equation 



V; 



f s (a 



>£<«) 



e£ T = 5gf 



(a) 



/.("<*) 



(4 Q) ) ; 



(13) 



(6) where due to the detailed balance of spin-flip scattering 

, does not depend on 



spin. 



From Eqs. (12 I and (13 I 



Co 



1 - P 2 ) Psfis + 1 



Co 

t(T) 



r 



0, (14) 



o-CD) /o-, els' 



S^)/S,sad Ai s = ^ (T) - 



where cr = cr(T) + crw, P 
C r(T) 1 S(T)+crW5Cl) ) P s = (S™ 
H^'. For metals the dimensionless parameter S 2 /Cq <C 1. 
By letting S 2 /Cq — ► the temperature gradient becomes 
constant and the conventional diffusion equations for the 
spin fj, s and charge accumulation fi c = (fifi^ + /i ) /2 are 
recovered ITOl: 



V 2 /i s (r) - ^ s (r) /l 2 si 



V 2 [Ac (r) + P/i, (r) /2] = 0, 



(15) 
(16) 



2 



withw = ((4 T) )" 2 + (4 i) )") 



-2\ -V2 



Let us consider a conducting magnetic wire with length 
A in the x direction. When we fix the charge and spin 
distribution functions at the left(right) ends to be h c l(r) 
and Msi(fl)) respectively, 



1 P 

McW = 2 (^ci + Mcfl) + J (fJ>sL + HsB.) 

P 



A 



(17) 



sinh (A/2 - x/l Bf ) 



sinh A 



Ms/? 



sinh (A/2 + x/l af ) 



sinh A 



(18) 



where A = A/Z s f measures the spin- flip scattering. When 
A 3> 1 and | / i B f — A| » 1 we find that fi s (x) = 



fi S Re x ^ ls! ) / sinh A is appreciable only in 



proximity of the edges. The vanishing of fx s in the bulk is 
independent of the thermoelectric forces applied by a tem- 
perature or chemical potential difference over the wire. By 
iteratively reintroducing the heat generation by spin-flip 
processes that are of order S 2 / Co we observe that then- 
effect is confined to the edges where V^a is significant. 
Gradients of the exchange potential and external Zccman 
magnetic fields do not appear explicitly at all. We con- 
clude that the spin accumulation cannot persist over dis- 
tances longer than the spin-flip diffusion length. Metallic 
spin diffusion therefore cannot explain the observed spin- 
Seebeck signal in permalloy with L/l s { ~ 10 5 [9]. 



3. Thermoelectric transport and thermal spin cur- 
rents 
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Figure 1: Spatial variation of the thermally induced spin accu- 
mulation (changes sign) and spin current in a magnetic wire for 
A = A// S j- = 5 (solid line) and 10. The scaling factor of the 
spin accumulation is ePgSAT/ (1 + PPg) and of the spin current 
<t(P- P')SAT/A. 



The local Fermi-Dirac spin distribution functions are 
defined by spatially and spin-dependent chemical poten- 
tials. Temperatures depend on position, but are not spin- 
dependent in the presence of the sufficiently strong inelas- 
tic scattering assumed here. In this Section we consider 
transport of the conduction electron spins in a magnetic 
metal experiencing thermoelectric forces, viz. an exter- 
nal electric field and gradients of the local chemical po- 
tentials and temperature. Using the Sommcrfcld expan- 
sion m Eqs. ((8J9J) we obtain for the local spin particle 

and heat currents Q°' = J^(r) — /i j' Q ^ (e, r), valid for 
C Q T 2 \d 2 a(e)\ £F \ «ct( £f ), 




rM (Q) /e 



5 T(i) 



S ( a ) 

£ T J V -V r T 

(19) 

conductivities obey the 
a C Ta^ in the limit 
and the t otal thermal conductivity k — 
= CqTct. Eq. (19 1 can be rewritten in terms 



The spin-dependent thermal 
Wiedemann-Franz law 



« 

A = 



of the charge, spin and heat currents, J c ( s ) = J^) ± J^' 



and Q = Q 



(T) 



Q ^ , respectively, 









■■( 








(20) 

in which P and P' stand for the spin-polarization of con- 
ductivity (c^Ief) and its energy derivative (d £ a^\ EF ). 

In the following we concentrate on the expressions for 
charge and spin chemical potential distributions induced 
in a one-dimensional magnetic wire by a constant temper- 
ature gradient, AT/A, which is part of an open electric 
circuit, i.e. in the absence of a charge current. The three 
unknowns /j, S L(R)&nd Afi are determined by boundary con- 
ditions. 

The chemical potential and temperature differences over 
the wire are defined as A/2 = [Ir — \Ili = Tr — T^. The 
open electric circuit boundary condition is J c (x) = 0. We 
consider two additional limiting boundary conditions: (i) 



3 



the reservoirs at the two ends of the wire are efficient spin 
sinks so that ^ s l(r) — 0; (ii) spin are not dissipated at the 
ends of the wire whatsoever, i.e. J s (— A/2) = J s (A/2) = 0. 
The reality might be somewhere between (i) and (ii). The 
heat current Q is uniform to leading order in S 2 /Cq. 

In the presence of ideal boundary spin sinks (i) the 
spin accumulation, Eq. (18 1, vanishes everywhere in the 
sample. The total chemical potential, Eq. (17), changes 
linearly in the sample and we find Ajl = eSAT, which 
implies a uniform thermally excited spin current J s = 
(P' - P)<rSAT/k when P' ^ P or equivalently P s ^ 0. 
This situation is similar to the classical thermocouple in 
which the wires of different materials (Sa ^ Sb) are elec- 
trically short-circuited at the ends. 

Without boundary spin sinks, (ii) J s (0) = J S (A) = 0, 
the spin accumulation distribution in the wire, Eq. ( |18[ ), 
becomes 



eP s SAT sinh(x/l sf ) 
(l + PP s )\ cosh (A/2) 



(21) 



changes sign when crossing the center, fi s (0) = 0, and its 
modulus increases exponentially at the edges (see Fig. [lj. 
At the ends x — ±A/2 we find finite spin accumulations 

eP s SAT , A 
/^ = -^ = - (1 + pPs)A tanh-. (22) 

Also this has an analogue with the classical thermopower 
effect [5] : a thermocouple consisting of wires from two dif- 
ferent materials, creates a voltage difference in the absence 
of currents (Ja(b) = 0). The spin accumulation has an ef- 
fect of the measurable thermopower S c = Ajl/ (eAT) = 
(S + ii S R/e) AT. The spin-flip scattering suppresses spin 
accumulation, but generates a thermally induced spin cur- 
rent 



_a(P-P') SAT ( sinh (x/l sf ) ' 
8 ~ A I 1 ~ cosh (A/2) 



(23) 



analogous to the electric current in a short-circuited ther- 
mocouple. When A> 1 the spin current equals that in (i) 
except at the edges, where it exponentially decays. In the 
opposite limit, in the absence of spin-flip scattering, the 
spin accumulation adopts a linear relation between con- 
tacts n s (x) cx xAT/A and the spin current vanishes. 

For efficient spin sinks at the edges (i), the spin cur- 
rent is uniform but the spin accumulation vanishes. The 
thermally induced spin accumulation in case (ii) is also 
suppressed sufficiently far away from the terminals. The 
spin accumulation is proportional to the temperature gra- 
dient and does not depend on length of the wire. A spin 
Seebeck effect should be observable in ferromagnetic wires 
in proximity of the edges. However, our results do not ex- 
plain the (approximately) linear inverse spin Hall voltage 
profile measured by Uchida et al. [5] in a permalloy thin 
film. 



4. Temperature dependence of chemical potential 
and spin density 

In this Section we discuss the temperature dependence 
of the chemical potential and spin density. The strong 
screening in metals enforces local charge neutrality inde- 
pendent of temperature. The associated variation of the 
chemical potential vs. temperature is sensitive to the en- 
ergy dependence of the density of states. 

The local electron density is the integral of the elec- 
tron distribution function over all cigenstates, for a nor- 
mal metal at thermal equilibrium n = /k (where 
/ = P = f*). Taking into account the energy depen- 
dence of the density of states, Af(e), the electron density 
in the limit (kgT/ /to) 2 <C 1. in which we may use the 
Sommerfeld approximation, reads |llj 



n= / deAT(e)fo(e;tio,T)=Ar(no)[lM) 



where 5n(0) = J Mo deAf(e) is a constant zero-temperature 
correction. SV = — e£ Td E lnjV(e)|^ can be interpreted 
as an intrinsic thermopower that does not depend on im- 
purity scattering. Recalling that a — e 2 MT> we observe a 
relation with the diffuse thermopower as S — SV + S-p. 
At thermal equilibrium the charge neutrality, n(T) — n(0) 
implies /i(T) = /io + eSj^T/2 , i.e., a temperature depen- 
dent chemical potential. The Maxwell identity dS/dn\T = 
—d\ijdT\ n = —eSj^/2 relates Sjy to the entropy <S. 

The above discussion can be extended to a ferromag- 
netic metal with spin-polarized density of states. The 
spin-dependent electron densities vary with temperature 
as 



< a >(M,T)= / deN^(e)h(e-^T) 



(25) 



Mo " + Sn^ (0) (26) 



Similarly, Sn^(0) = J" deAf^(e) 7 where 5 



(a) 
'AT 



-eA/T<9 e lnA/' (Q) (£)| W) and ^ {a) is the density of states 
at the chemical potential. The transport thermopower 
is S^ — S& + . At finite temperature, the spin- 
polarized eSj^ T/2 act as intrinsic thermoelectric poten- 
tials, resembling the Zeeman splitting of energy bands by 
an external magnetic field. The charge neutrality condi- 
tion in the normal metal holds in the fc rrom agnet as well, 
i.e. d(n^ + nW)/dT = 0. Using Eq. ^ this condition 
implies 



ia Mo _Arms$ 



dT 



(27) 



which reduces to that in the normal metal when J\f(^ — 
J\f^K A change in the equilibrium temperature thus shifts 
the chemical potential (the Fermi level) just as a gate volt- 
age but with magnitude proportional to the intrinsic ther- 
mopower ((5/i = eS^fST). A similar effect exists in the 
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presence of a uniform magnetic field, where the charge neu- 
trality condition implies a magnetic field-dependent chem- 
ical potential in the fcrromagnet known as the magneto- 



Coulomb effect [Hrrg. Substituting Eq. ((27) into the ex- 
pressions for the spin-polarized electron densities, we find 
for the variation of the spin density n s = — as a 
function of temperature as dn s /dT = eNSj^{P^f — Pj^) 
or, after integration over temperature, 



n.(T) 



n.{0) = ^NS V T{P V 



K) 



(28) 



in which P^ and Pj^ are the spin polarization of the den- 
sities of states (J\f( a )) and its energy derivative (d^Af), 
respectively, and Af = AfW + Af^>. The spin density is 
thus temperature dependent when Ps^ = (Pjj — Pj/) /(l — 

The density of states in a ferromagnet is in general also 
temperature dependent, dAf/dT = (dAf/dA) (dA/dT) ^ 
0. At thermal equilibrium the spin density and the ex- 
change splitting are closely related to each other A(T) cx 
n s (T). In the Appendix we show that even far below the 
Curie temperature this temperature dependence may be 
significant. 



5. Thermally induced spin density/polarization 



In Section 3, we derived the spin accumulation in a fer- 



romagnetic wire under a temperature gradient, Eq. (21) 



The chemical potential is a function of the electron density 
and the temperature, as shown in the previous Section. A 
local variation of the spin chemical potential from equilib- 
rium 5/j,( a '(x) = fj,( a '(x) — [Xq can be expanded in terms of 
deviations of non-equilibrium spin-polarized electron den- 
sities and temperature (rS a \ T) from their equilibrium val- 

ues (n a) ,T ) 



5^ a \x) = 



8n {a \x) 



V dT 



ST(x). 
(29) 



The derivatives can be found by an extension of Eq. ( 26 1 



by replacing /o with the local spin-dependent distribution 
functions f a \e;^ a \T), as 



dnF) 
dT 



1 



eS 



(«) 



(30) 
(31) 



In a bulk ferromagnet with spin-flip scattering we have 
seen that /j, s (x) — 0, such that 



5n m (x)/Af m -5n a) (x)/Af a) 



e(S$P-S%>){T(x)-T ), (33) 



which is a spin density difference (spin polarization) in- 
duced by the temperature gradient. Eq. ((33), regardless 
of its position dependence, is equivalent to Eq. ((28) for 
the temperature dependence of the spin density at ther- 
mal equilibrium. The spin polarization, Eq. (33 1, adopts 



a linear profile and changes sign in the center of the mag- 
netic wire, where To = (Tl + Tr)/2. It originates from the 
temperature dependence of the chemical potential in the 
presence of spin polarization of the intrinsic thermopower 
in the ferromagnet. 

The temperature dependent spin polarization and lo- 
cal exchange splitting in the bulk ferromagnet is a local 
equilibrium effect that cannot drive a non-equilibrium spin 
current into lateral contacts. The observed spin signal in 
the experiments by Uchida et al. [9] can therefore not be 
due to the spin polarization induced by the intrinsic ther- 
mopower. A thermally induced spin pumping mechanism 
due to magnon propagation in the bulk ferromagnet [14] 
may explain the observed spin signals. 

6. Conclusions 

We presented a theoretical analysis of thermoelectric 
transport in ferromagnetic metals using an extension of the 
spin-dependent Boltzmann and diffusion equations. We 
showed that in a bulk ferromagnetic metal subject to a 
temperature gradient a pure thermal spin current is gener- 
ated by spin-flip scattering whereas the spin accumulation 
vanishes. The temperature dependent spin density aris- 
ing from an energy-dependent density of states generates 
a spatially varying local spin polarization along the ferro- 
magnetic wire. The spatial dependence of the exchange 
splitting induced by a temperature gradient to leading or- 
der does not lead to effects on spin transport. This work 
could help to better understand thermoelectric phenomena 
in metallic ferromagnets. 
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In a ferromagnet subject to a temperature gradient, the 
spin chemical potential thus contains two contributions 

H„{x) = 5n^(x)/Af m - 6n^(x)/Af^ 

+ e(S { ^-S^)(T(x)-T ) (32) 



Appendix: Temperature dependence of the ex- 
change potential in the Stoner model 

Here we analyze the temperature dependence of the 
chemical potential and equilibrium spin densities based on 
the Stoner theory of ferromagnetism. The Stoner criterion 
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W(/io) J = 1, in which J is the Stoner exchange parameter, 
describes the onset of ferromagnetism. A stable magnetic 
ordering leads to an exchange splitting of energy bands 
for opposite spin directions, which is related to the spin 
density by A(T) = Jn s (T). J is an essentially atomic-like 
exchange integral that does not depend on temperature. 
Here we study the temperature dependence of the chemical 
potential and the spin density self-consistently. 

We emphasize that the Stoner model is woefully inad- 
equate at elevated temperatures since it completely ne- 
glects spin wave excitations of the magnetic order parame- 
ter. Consequently, the critical temperatures are drastically 
overestimated by the Stoner criterion. Nevertheless we 
believe that at temperatures sufficiently below the Curie 
transition, the Stoner model, as a simple implementation 
of density-functional theory, can provide useful qualitative 
insights. 

The density of states depends on temperature via the 
exchange splitting. In the rigid-band model Af( a > = 
Af(e + aA). The charge neutrality condition d(n^ + 
n- 1 ) /dT = now implies 



, , 3M Kr df 







(A-1) 



where Af = Af' +AA and / = /I = at thermal equilib- 
rium. Using the identities 

df _ ( df\ ( e - n ( 8(jl\ , ( df 
dT 



Be 



T 



dT 



dT 



(A-2) 



and dAf/dT = (dA/dT) (dAf/dA) , we find for the tem- 
perature dependence of the chemical potential 

d V -1 f , ( e ~ aA K r, m\ ( d f" 



dT 



-x 



de 



e — /i 
T 



■A/"(e,T) 



de 



X 



df 



dN(e,T) 
de—g^m 



(A-3) 



where x(T) = / deAf(e,T) (-d e f) is the (Pauli) suscepti- 
bility. The first term above is the intrinsic thermopower 
(= eStf), and the second term is the correction by the 
temperature dependence of the exchange splitting. The 
above equation can be rewritten as 



dfi 

df 



X^dA 
X dT 



(A-4) 



in which Xs — X^ ^ X 1 ■ 

Next we are interested in the variation of spin density 
by temperature. The self-consistency condition reads 



dn s {T) 
dT 



1 dA(T) 
J dT 



(A-5) 



With Eqs. ( A-2|A-4 ) , the temperature dependence of the 
spin density becomes 



dn s (T) 
dT 



1 - 



e(5V - S^Jxs 
Jx (i - (Xs/xf) 



(A-6) 



in which Sj\f s and \s are obtained as intrinsic thermopower 
and susceptibility but with substituting Af(e,T) — > 
Af s (e,T). Equation (A-6 1 is exact within the limits of the 
validity of the rigid-band Stoner model. Neglect of the in- 
teraction correction dn s /dT sw e(SV — Sj\f,)Xs is equivalent 
to our former result Eq. ( |28[ ). At the Curie temperature 
J\ = 1 and dn s (T)/dT diverges. At low temperatures, 
assuming that J\ does not change much from the Stoner 
criterion 



dn s (T) e(P M - P' u ) X Sm 
dT ~* l-Jx(l-P^) 
^ (P M - P' M ) 



p2 



(-eC T) d £ Af\, 



(A-7) 
(A-8) 



The temperature dependence of the gap follows from Eqs. 
( A-6|A-5 1 , which leads to contributions to the transport 
thermopower as explained below Eq. (19 1. 
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